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We analytically derive a more accurate formula for the power spectrum of the biased objects with 
the primordial non-Gaussianity parameterized not only by the non-linearity parameter /nl, but 
also by <7nl and tnl which characterize the trispectrum of the primordial curvature perturbations. 
We adopt the integrated perturbation theory which was constructed in Matsubara (2011) We 
discuss an inequality between /nl and tnl in the context of the scale-dependent bias, by introducing 
a stochasticity parameter. We also mention higher order loop corrections into the scale-dependency 
of the bias parameter. 
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Primordial non-Gaussianitics have been widely well-studied as new probes of the mechanism of generating primordial 
curvature perturbation, e.g., inflationary physics. The observation of the cosmic microwave background (CMB) 
(-h \ anisotropics is one of the best ways to hunt for the primordial non-Gaussianities through the measurement of the bi- 
Q_j and tri-spectra of the temperature anisotropies. For example, a current best limit for the primordial non-Gaussianity 
is obtained as — 10 < /nl < 74 at 95% confidence level (CL) [H with /nl being a local-type non-linearity parameter, 
which characterizes the bispectrum of the primordial curvature perturbations. The standard single field inflation 
model predicts that /nl is the order of 10 -2 and hence if we would detect that /nl > 0(1) it implies that our 
. universe has not experienced the standard inflationary stage with a single scalar field. 

Thanks to the progress of the cosmological observations, recently, the large scale structure (LSS) of the universe 
will also bring us the fruitful information about the primordial fluctuations. While in the observation of LSS the 
understanding of the gravitational non-linear evolution of the density fluctuations is more important on smaller scales, 
| the primordial non-Gaussianities affect the clustering of the biased objects on large scales. It has been found that the 
Q\ . local- type primordial non-Gaussianity, /nl, produces significant scale-dependence on the bias parameter, Ab ex 1/fc 2 . 
The constraints on the primordial non-Gaussianities obtained from the observations of such scale-dependence of the 
bias are comparable to those from the CMB observations, e.g., —31 < /nl < 70 at 95% CL in Ref. [3(. 

However, in order to make ready for the upcoming precise observational cosmology, on the theoretical side, more 
deeply understanding of the effect of the primordial non-Gaussianities should be necessary. There are a lot of literatures 
about the precise formula for the scale-dependent bias both analytically and numerically and also the analysis of the 
higher order effects. Recently, Matsubara (2012) Q has shown an accurate formula of the scale- dependent bias with 
primordial non-Gaussianity parametrized by /nl, by making use of the integrated Perturbation Theory (iPT) [l[. 
By adopting this formula, it is not necessary to use the high peak approximation and also the peak-background split 
picture, both of which have been considered as useful tools to derive the bias parameter in the case with the primordial 
■ non-Gaussianity. 

Together with deriving the precise formula of the bias parameter for /nl case (see, e.g., it should be important 

to investigate the effects of the higher order primordial non-Gaussianities, that is, the trispectrum of the primordial 
curvature perturbations. The non-linearity parameters denoted by <?nl and tnl were introduced to characterize the 
amplitude of the primordial trispectrum. There are also lots of literatures which study these non-linearity parameters 
theoretically and observationally. One of the interesting theoretical issue about these higher order non-Gaussianitics 
is an inequality between /nl and tnl, so-called, Suyama-Yamaguchi inequality 0-0 ■ If the equality tnl = 36/^/25 
is justified, the primordial curvature perturbations should be sourced from a single scalar field. The case with 
Tnl > 36/^ L /25 indicates that the primordial curvature perturbations might be sourced from multi-scalar fields. 
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There have been several works about the effects of these higher order parameters c/nl and tnl on the scale- 
dependency of the bias parameter. In Ref. (Ioj |. the authors derive a scale-dependent bias in the case with non-zero 
3nl and find that <7nl gives the same scale-dependency as the case with /nl- They also give a constraint for 
by using the same data set as in Ref. Q. For tnl, although there are not any available constraints, Refs. fill Il2| 
analytically derive the formulae of the scale-dependent bias with tnl (Refs. (HEI 

as related works). Numerical 

analysis for the effects of such higher order non-Gaussianities on LSS has been shown in Refs. | lSj - flSl ] . However, 
comparing with /nl case, the precise formulae for <?nl and tnl cases might be less understood. 

In this paper, we extend the analysis of the scale-dependent bias in Ref. [|| to the case with higher order non- 
Gaussianities characterized by the non-linearity parameters tnl and <?nl, by making use of iPT. In terms of iPT, the 
effects of <7nl and tnl appear in two-loop order contributions and hence we also investigate the other contributions 
from the higher order loops. 

This paper is organized as follows. In the next section, we derive a power spectrum of the biased objects with 
including the higher order primordial non-Gaussianities, by making use of iPT. In section InTl we focus on the leading 
order contributions of ^nl and tnl, and discuss the scale- and redshift-dependence of the bias parameters. We 
introduce the stochasticity parameter in order to investigate the inequality between /nl and tnl- In section ITVl we 
discuss the higher order effects which are linearly proportional to the non-linearity parameters. Section [V] is devoted 
to summary and discussions. We plot the figures of this paper with adopting the best fit cosmological parameters 
taken from WMAP 7-year data @. 

II. POWER SPECTRUM OF THE BIASED OBJECTS WITH THE PRIMORDIAL 

NON-GAUSSIANITIES 

Here, we focus on the local-type primordial non-Gaussianity because the several literatures have shown that other 
types of non-Gaussianities such as equilateral-, orthogonal- and folded-types are not so effective on the structure 
formation of the Universe, compared with the local type. 

A. Local-type non-Gaussianity 

Let us introduce the local-type non-linearity parameters. The primordial non-Gaussianities of the curvature per- 
turbations, $, are usually characterized by the higher order spectra as 

($(k)$(k')> = (2 7 r) 3 ,5 (3) (k + k') J P*(fc), 

(*(ki)$(k a )*(k3)> = (2^) 3 ,5 (3) (k 1 + k 2 + k3)B*(ki,ka, k 3 ), 

($(k 1 )$(k 2 )$(k 3 )$(k 4 )) c = (2^ 3 ( 5( 3 )(k 1 +k 2 +k 3 +k 4 )T $ (k 1 ,k 2 ,k3,k 4 ), (1) 

where a index, c, denotes a connected part of the correlation functions. For the local-type case, the higher order 
spectra are parameterized as 

Tb(/ci,fc 2 ,fc 3 ,fc 4 ) = 63NL [P^(ki)P^(k 2 )P^(k 3 ) + 3 perms.] 
25 

+ yr NL [P <& (fc 1 )P $ (fc 2 )P$(|k 1 + kg|) + 11 perms.] . (2) 

In the above expression, we call scale-independent parameters, /nl, 3nl and tnl, as the local-type non-linearity 
parameters. As the simplest local-type model, we consider that the primordial fluctuations are expanded as 

$(x) = $ g + /nl(*g(x) 2 -($ g (x) 2 ))+. 9 nl$g(x) 3 + -- - , (3) 

where <E>g denotes pure Gaussian fluctuations. We have a special inequality for the local-type non-Gaussianity which 
is given by 0-0] 

TNL > 7rr/NL- (4) 
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For the case where the primordial curvature fluctuations are sourced from the quantum fluctuations of a single scalar 
field and is given by Eq. we can realize the equality of the above relation. If there exist the multiple scalar 
fields in the early Universe, which are seeds of the primordial curvature fluctuations, we have tnl > 36/25/^ L . For 
example, in the case with 

$ = &j(x) + V; G (x) + /* L (0 G (x) 2 - (0 G (x) 2 )) , (5) 

we have 

fL _36 



/NL 



and 



(1+0 2 ' 25(l + 3 

TNL = + 
3fi 

> ^/nl, (for £ > 0) (7) 

where we have used = and £ = P^/P^ with P^(i = </>, ?/>) being the power spectrum of each component. Hence, 
because of this fact, it is important to investigate the relation between tnl and /nl for studying the origin of the 
primordial fluctuations. 



B. The power spectrum of the biased objects in integrated perturbation theory 

Let us derive a formula of the power spectrum of the biased objects with the primordial non-Gaussianity by making 

use of iPT. Following Refs. jl|, 0] , we introduce multipoint propagators T£ of the biased objects X in "Eulerian" 
space which is given by 

/ ^"<M k ) \ - f^-\3-3n X (3)f^ j i \r _l^r(") 



^ L (ki)---<W L (k„ 



(27r )3-3n j( 3) (ki + ... +kn _ k)r («) (ki; ... jkn)5 (g) 



where ^L(k) and Sx(k) are respectively the Fourier transforms of the linear density field and the number density field 
of the biased objects in "Eulerian" space. The linear density field ^(k) is given by 

<S L (k) = M(fc)$(k). (9) 

Here, the proportional factor Ai(k) is given by 

M(k) _ ^ D(z) k*T(k) 

where T(k), D{z), H$ and O m o are respectively the transfer function, the linear growth factor, the Hubble's constant 
and the matter density parameter. Here z* denotes an arbitrary redshift in the matter-dominated era. 
The power-, bi- and tri- spectra of ^(k) are respectively given by 

P L (fc) = X(fc) 2 P$(fc), 
S L (ki,k 2) k 3 ) = ^(fc 1 )M(fc 2 )X(fc 3 )P 4 ,(k 1 ,k 2 ,k 3 ), 

^(k^ka.ks,^) = ^(fc 1 )M(fc 2 )>l(fc3)A4(A ; 4)T $ (k 1) k 2) k3 J k 4 ). (11) 

Up to the leading order of the non-linearity parameters /nl, 3nl and tnl, in terms of the multipoint propagators, 
the power spectrum of the biased objects is given by 

Px(fc) = [rW(k)] 2 p L (fc) + r^(k) J ^rf (p,k-p)B L (k,- P) -k + p) 

+ir^(k) J d ff^ rff (pi, p 2 , k - pi - P2 )r L (k, - Pl , - P2 , -k + Pl + P2 ) 
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d 3 p 1 d 3 p2 r ,(2) l 



r^( P i,k - P i)r^(-P2, -k + P 2)T L ( Pll k- Pl , -p 2 , 



(27r) e 



P2j 



(12) 



While the first line in the above equation has been already obtained and investigated in details in Ref. Q, here, we 
look more deeply into the second and third lines. 

The multipoint propagators include the non-linear evolution of the matter density field and also the bias functions. 
In order to derive a more useful expression for the power spectrum, let us introduce renormalized bias function in 
" Lagrangian" space following Ref. |4| , which is defined as 



C L(k 1 ,k 2 ,...,k n )^(2 7 r) 3 "| ^ 



5 n S L x (p) 



(13) 



where 8\ is the number density field of the biased objects in "Lagrangian" space. The multipoint propagators of the 
biased objects on large scales (ki + k 2 + ■ • • + k„ — > 0), where the non- linear evolution of the matter density field is 
negligible, are written in terms of c 1 ^ as 



r£(k) l + c^k), 



4 n) (k 1; k 2 ,. 



> kfi) 



c^(ki,k 2 , • • • ,k„) for n > 2. 



(14) 



From Eq. (fTB"]) . the renormalized bias function is calculated once the relation between the linear density field and the 
number density field of the biased objects in "Lagrangian" space is determined. Based on the Press-Schechter (PS) 
picture, the mass function is determined by the probability that the value of the "smoothed" linear density field over 
mass scale M exceeds a critical value S c , denoted by P(M,S C ). The smoothed linear density field over mass scale M 
is given by 



d 3 k 

(27p 



W(kR)6 h (k), 



(15) 



where R is a comoving scale corresponding to the mass scale M and W(kR) is a usual top-hat window function. 
Following Ref. [|J, the renormalized bias function can be written as 



C n(k-l7 'j k„) — 



(-l) n 7& l ^&'^ WfaR) ■ ■ ■ W(k n R) 



dP(M, S c )/dM 
Assuming a universal mass function, wc have 

dP(M,S c ) f M F(v)d]na M 



(16) 



dM 
d n P(M,S c ) 



2 dM ' 
(-!)>- l)!^ 1 (-1)- {m) 



2(5" 



(17) 



rn=0 



where /mf(^) with v = 5 c /ctm is a multiplicity function. Then, using the above expressions, the renormalized bias 
functions are reduced to 



c£(fci,-- M) 



(n-l)! 



E 

\m=0 



(-1)* 



^/mfV) )w(hR)---W(k n R) 



b^{M)W{k 1 R)---W{k 2 R) 



(n-l)! 



din ctm 



[W(kiR)---W(k 2 R)] 



(18) 



where the scale independent function, b^(M), given by 



b-(M) = ( — 



&M J /mf(i / ) 



(19) 
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We can take into account the effect of the primordial non-Gaussianity on the renormalized bias function by using 
the non-Gaussian multiplicity function. We discuss the non-Gaussian corrections to the renormalized bias function 
in a later section. For example, wc have 



/mf(^) = /psO) = 




in the original PS theory, and 

J-qve-^l\ (21) 
for the Sheth-Tormcn (ST) fitting formula with p = 0.3, q = 0.707 and A(p) = [1 + TT-^ 2 2- p T(l/2 - p)]" 1 . 



/mf(i>) = /stM = A(p)\l~ 



iqv 2 )P 



III. SCALE-DEPENDENT BIAS WITH PRIMORDIAL NON-GAUSSIANITIES 



Here, by using the simple expressions for the multipoint propagators on large scales, which are given in the previous 
section, let us write down the bias parameter with respect to the non-linearity parameters /nl, <7nl and tnl- Defining 
a bias parameter as 



P x (k) = b 2 x (k)P h (k), 



(22) 



on large scales (k — ¥ 0), we have 



b 2 x (k) a b x {Mf + Ab 2 x (k) 



Ab x (k) = 4/, 



XL 



h(M) f d 3 Pl 



h(M) fd 3 Pl d 3 P2 



M(k) J (2tt) 3 

DffNL + TT'TNL 

y 

XCg(-pi, -p 2 ,Pl + P2)M(j>l)M(p2)M{\pi +P2|)P*(P1)P*(P2) 



M(k) J (2tt) 6 



25 1 



d 3 p l d 3 p 2 l , \ L / , D/ , D , , 

— J^yT c 2(Pi,-Pi) c 2(P2,-P2)^l(pi)Pi> 2 )- 



(23) 



where b\(M) = 1 + c L is the scale- independent Eulerian linear bias parameter and Ab x (k) denotes a scalc-dcpcndcnt 
part due to the primordial non-Gaussianity. The renormalized bias functions up to the third order are given by 



C2(-Pl,Pl) 



s 2 



W( Pl R) 



dW( Pl R) 
din ctm 



c£(-pi,-pa,Pi + Pa) = b%(M)W(piR)W(p2R)W(\pi+p2\R) 



2 + 26 c (6i(M) - 1) + 5 2 c b\{M) 



Substituting these expressions into Eq 

bi(M) 



Ab x (k) » 4/, 



NL 



M{k) 



%(M) + 2 



we have 

l + <f c (&i(M)-l) 



diner m 



[W(piR)W(p2R)W(\pi + p 2 \R)} ■ (24) 



S 2 



M 



25 



,9nl + 27 TNL 



25 



i 



9 TNL M{k) 2 



MM) 
M(k) 

2 + 25 c (b 1 (M)-l) + 8 2 c b\{M) 
1 + S c (h(M) - 1) 



i | d\nS 3 (M) \ 
In a M J 



*mS 3 (M) 



b\{M) + 2 



S 2 



(25) 
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where <jm is a variance of the smoothed density field, §m, and S^M) is a skewness given by 



S 3 (M) 



6 



d 3 p 1 d 3 p 2 
(2tt) 6 



^(pxiijWCpaiZjWdpi + P2\R)M( P i)M(p 2 )M{\p 1 + p 2 |)P*(pi)P*(P2). 



(26) 



The first line in the above expression has been already shown in Re f. Lil and gives the scale-dependence with Ab(k) oc 
1/k 2 . A recent paper by Biagetti, Desjacques and Riotto (20f2) [1J] have shown a non-Gaussian correction to the 
bias including the effect of tnl and we find that their expression of the non-Gaussian correction coincides with the 
second line in our Eq. (|25|) . by employing the PS mass function. Gong and Yokoyama (2011) [ll| derived the similar 
correction term dependent on tnl as the last term in Eq. (|25|) by making use of the high peak limit formalism. By 
employing PS mass function and taking high peak limit, we find that the last term in our expression would agree with 
the TNL-correction term in Ref. fllT ] . In this sense that the results in the previous works are included in Eq. (|25l) as 
limiting expressions, this expression is a general expression of the scale-dependent bias depending on the non-linearity 
parameters, / NL , #nl and r NL . 



/nl vs #nl 



Let us focus on gNL-corrcction term. From Eq. (j25[) . we find the c/NL-correction gives the same scale-dependence 
as the /NL-correction. In this sense, the non-linearity parameters c/nl and /nl are degenerate with each other in 
the observation of the biased objects. One way to distinguish these two non-linearity parameters is to see the mass- 
dependence (6i (M)-dependence) of the non-Gaussian corrections, which has been shown in Ref. [l8|. Here, we focus 
on the redshift-dependence at fixed scale of the non-Gaussian corrections. In Fig. [TJ we plot Ab as functions of the 
redshift with fixing k = 0.01 h Mpc -1 and M = 5 x 10 13 /i _1 M . We set /nl = 40 for red circles and g N L = 5 x 10 5 for 
blue boxes. From this figure, we find that the redshift dependence of the scale-dependent part of the bias parameter 
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FIG. 1: Ab 2 x as functions of the redshift with fixing k = 0.005 h Mpc" 1 and M = 5 x 10 13 ft _1 M . We set / NL = 40 for red 
circles and jnl = 5 x 10 5 for blue boxes. 

is different between the cases with non-zero /nl and (?nLj and also we can see that the higher redshift observations 
would be a powerful tool to obtain a constraint for (?nl- 



B. 



/NL VS TKL 



Next, let us focus on the TNL-correction terms. As we have mentioned, a significant constraint on the tnl parameter 
is a powerful probe of the origins of the primordial adiabatic fluctuations, since it has a special inequality with /nl- 
To focus on the inequality between /nl and tnl in the context of the scale-dependent bias, we introduce the power 
spectrum of the matter density field, P m , and also the cross power spectrum, P m x- On large scales where the 



non-linear evolution of the matter density field is negligible, these power spectra are approximately written as 

Pm(k) « P L (fc), 

P L (fc) f d 3 Pl 



, 25 ^ P L (fc) / d 3 pirf 3 P2 



x4 , (-pi,-p 3 ,pi+p a )M(pi)^t(pa)M(|pi + P3|)P*(pi)P#(p3). (27) 



9 "V A4(fc) 7 (2tt) 6 

By making use of the matter power spectrum and the cross power spectrum, we estimate a stochasticity (see, e.g., 

2 



<GO= - ( --~o-«£>= )' 



FIG. 2: The diagrammatic picture for the stochasticity parameter f(k) — 1. The double solid line, the dashed line, the solid 
line and the circle represent respectively 5x , 5m, 5l and the multipoint propagator. 




FIG. 3: The stochasticity parameter as a function of k with fixing z = f , M = 5 x fO 13 h 1 Mq and /nl = 40. The red thick 
line is for the case with tnl = 2 x 36/nl/25. The red dashed line is for tnl = 5 x 36/nl/25 and the red dotted-dashed line for 
r NL = 10 X 36/£ L /25. 



[E EH, H3] ) ■ Here, for convenience, we define a stochasticity parameter as 

_ P m (k)Px(k) 



Pmx(fc) 2 ' 



At the leading order, we can obtain 



25 



r(k) ~ 1 + ( -r NL - 4/^ L 



6i(fc) 2 A4(fc) 2 



° 4(p,-p)P L (p) 



(2tt 



(28) 



(29) 



1 In Refs. [ll I lSi, [20t . the stochasticity parameter is defined as r(k) = P m x(k) / \J Pm(k)Px (k). In Ref. Il5l| . the authors introduce 
r(k) = P x (k)/Pm{k) - (P m x(k)/Pm(k)) 2 and Ref. [JJ] gives x(fc) = PmxMVC^xW-fmW)- Note that these parameters have the 
similar dependence on the non-linearity parameters although each definition is slightly different. 
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The diagrammatic picture of the stochasticity parameter f(k) — 1 is shown in Fig. [2] following the rules given in Ref. 

In this figure, the double solid line, the dashed line, the solid line and the circle represent respectively Sx , S m , <5 L 
and the multipoint propagator. In Fig. we plot the stochasticity parameter as a function of the wavenumber, k, 
with fixing z = 1, M = 5 X 10 13 /i _1 M and / NL = 40. The red thick line is for the case with r NL = 2 x 36/ NL /25. The 
red dashed line is for tnl = 5 x 36/^/25 and the red dotted-dashed line for tnl = 10 X 36/^ L /25. From this figure, 
we find that due to the deviation from the equality as tnl = 36/^ L /25 the stochasticity parameter f{k) deviates 
from unity on large scales. This means that the observations of the stochasticity parameter on large scales could be a 
powerful tool to study the inequality between /nl and tnl • This fact was also discussed in Ref. [13l Il5j for two-field 
inflation model with different parameterization. Note that the stochasticity includes various uncertainties due to the 
higher order corrections [l[ and also astrophysical processes. The above discussion about the stochasticity given by 
Eq. ([2T>]) is ideal and hence the detection of the effect of r NL needs more careful investigation. In later section, we 
discuss the higher order corrections into the stochasticity. 



IV. HIGHER ORDER CORRECTIONS 



In the context of the iPT, the correction terms from the non-zero (?nl and tnl discussed in the previous section are 
two-loop order corrections, while the /NL-correction term is the one-loop order. Of course, we have other higher order 
corrections which we have neglected in the previous section. Although it is difficult to check all the contributions 
analytically, in the large scale limit where the non-linear evolution of the matter density field is negligible, the 
estimation for such contributions becomes somewhat easier. Here, as the higher order corrections, we focus on the 
corrections in the pure Gaussian fluctuations and also those which are linearly proportional to the non-linearity 
parameters and discuss in what case we do not need carefully to treat such higher order corrections. 



A. Non-Gaussian corrections in multipoint propagators 



As we have discussed, the multipoint propagators on large scales are rewritten in terms of the renormalized bias 
function as Eq. (|14[) and the renormalized bias functions are given by Eq. (fTSJ) in terms of the multiplicity function, 
/mf(^)- In case the primordial non-Gaussianity exists, it has been known that the multiplicity function is modified 
(see, e.g., Ref. [H|) and hence the multipoint propagators should be also dependent on the non-linearity parameters. 
To include the corrections due to the primordial non-Gaussianity, we introduce a correction function as i? and 
rewrite the multiplicity functions as /mf = /ps(st) x -R NG - There are several works about i? NG analytically and 
numerically. For example, based on the Edgcworth expansion formula, up to the terms which are linearly proportional 
to the non-linearity parameters, we have 



pNG 

-"■Ed 



f 



« 3 (M) 
6 



1 dn 3 {M)/dM 
v 69 In om I dM 



H 2 {u) 



( Kj(M) 

V 24 



H 4 (u) 



1 dm{M)/dM 
v2±d\\\<T M /dM 



H 3 {v) 



(30) 



with k„ = (ijj) c /aj f and H n (x) being Hermite polynominal. The higher order moments k 3 and K4 are respectively 
proportional to the non-linearity parameters /nl and Snl (also tnl)- As another example, Matarrese, Verde and 
Jimenez (2000) [2l| has shown a fitting formula for i?NL which is given by 



i? NG 



MVJ 



v k 3 (M) dK 3 (M)/dM 
fa 6 6d\na M /dM 



v 2 k 4 (M) dK4,{M)/dM 



12 249 In a M /dM 



exp 



k 3 (M) 3 , k 4 (M) 



G 



24 



(31) 



with 8 n = yj\ - 2v n -' 2 n n (M)/n\. 

Let us evaluate the amplitudes of the modifications of the renormalized bias functions due to the primordial non- 
Gaussianity. Hereinafter, we adopt the Shcth-Tormcn fitting formula as a Gaussian mass function and the Edgcworth 
expansion formula as a non-Gaussian correction of the mass function. In this case, the first order renormalized bias 
function on large scales (k — >• 0) is 



c i,ng(^) 



1 



qv 1 - 1 
Ks(Af) , 



om) /mfH 
1 + {qv 2 )P 

v _ K 4 (M) 

V ' 6 



-1 
cm 



/ps(") 
/ps(f) 



PNG 
■"-Ed 



o om<TM \ v 



1 8k 4 (M) 
12 d In a m 



(32) 
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In Fig. |4l we plot the ratio of the bias parameter in Eulcrian space b\ (M) = l + c\ with the primordial non-Gaussianity 
to that in the Gaussian case, as a function of v. The red circles are for the /nl = 40 case, the blue boxes for the 
5nl = 5 x 10 5 case and the green diamonds for the tnl = 10 x (6 x 40/5) 2 . From this figure, we find that in the 
range with 1.5 < v < 4.5 corresponding to the mass range of the biased objects with 10 12 < 7\//[/i _1 A/©] <5x 10 14 
for z = 1, the contribution of the modifications of the mass function by the primordial non-Gaussianity to the bias 
is negligibly small. For /nl = 40 and v = 4.5, the effect of the primordial non-Gaussianity is about 4%. For the 
higher peak objects (y > 5) the effect of the primordial non-Gaussianity is expected to become larger and hence for 
much higher peak objects we could not ignore the contribution any more. Such kind of non-Gaussian effects can be 
obscrvationally seen in the number counts of the very massive objects which are minor components in the observations 
of the clustering of the biased objects. 




1.5 2.0 2.5 3.0 3.5 4.0 



FIG. 4: The ratio between the non-Gaussian first order bias parameter and the Gaussian one, b x (M)/b% (M) as a function 
of v. The red circles are for the /nl = 40 case, the blue boxes for the $nl = 5 x 10 case and the green diamonds for the 
tnl = 10 x (6 x 40/5) 2 . 



B. Higher order loops in terms of the multipoint propagators 

In the previous subsection, we find that the modification of the multipoint propagators due to the primordial non- 
Gaussianity is small for not so high peak objects. Hereinafter, we consider the case with M — 5 x 10 13 h~ 1 M Q at 
z = 1.0 and neglect the effects of the primordial non-Gaussianity which appear in the multipoint propagators. 

Next, let us focus on the higher order loop corrections in the power spectrum of the biased objects, in terms of 
the multipoint propagators, . The power spectrum of the biased objects, up to the terms which are linearly 
proportional to the non-linearity parameters, is obtained by 



Px(k) = £ 



/ dV (2 ^ ) l P "" 1 4 l) (pi,P2,--- ,k-P n _ 1 ) 2 abl)PLb 2 )---PL(|k-P„_ 1 |) 
[)\ J (2tt) 3 " iy( P i,--- ,p n -i,k-P n -i) 



Xr£ +1) (- Pl , • • • , - P „, -k + P„)P L (pi) • • • ft(p„-l)S L (k - Pn-1, -Pn, -k + Pn) 



1 f d 3 px ■ ■ ■ d 3 p n+1 ( 
3(n- 1)! J (2tt) 3 »+ 3 



+ o;„ TYi / ,o,^n+3 r x (PV ,Pn-l.k-P„-i) 



xr^ +2) (-pi, • • ■ , -p„+l, -k + P„+i)P L (pi) • • ■ P L (p, l -l)T L (k - Pn-l, -Pn, -Pn+l, -k + P„+i) 
d 3 p! ■ ■ ■ d 3 p„-lrf 3 g-p(n) , 



n=2 



4(^2)! J (2^ ^ (P1 '"' -P«-l' k - P n-lJ 
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xT x (-Pi) ■ • • > -Pn-2, -q, -k + P„- 2 + q) 

xP L (pi) • • • P L (p„-2)T L (p„-i,k - P n _!, -q, -k + P„_ 2 + q), (33) 
where P„ = pi + • • • + p n . In the right hand side of the above equation, the first term in the bracket appears 




FIG. 6: Diagrams linearly proportional to the primordial bispectrum. 




FIG. 8: Other diagrams linearly proportional to the primordial trispectrum. 



even in the case with pure Gaussian fluctuations, corresponding to Fig. [5] The second term in the bracket represents 
the contributions which are linearly proportional to the primordial bispectrum parameterized by /nl- This term is 
corresponding to Fig. |B1 Both the third and last terms are the contributions which are linearly proportional to the 
primordial trispectrum parameterized by <?nl and tnl- These terms are respectively shown by Figs. 0and|H] Such 

(n) 

kind of higher order corrections for the power spectrum of the matter density fields arc estimated by changing T x to 
r^' in the above expression. As we have mentioned, the non-linearity of the matter density fields becomes negligible 
and r2' — > 1 and r„' — > (n > 2) on large scales (fc — > 0). Hence, for the power spectrum of the matter density 
fields such higher order correction terms are negligible on large scales. 

Using Eq. (|14[) . the power spectrum of the biased objects including the higher order corrections is reduced into 

Px(k) « b^kfP^ + AfMk)^^! |^(-Pi,Pi)ft(pi) 

+ (6 flNL + T 7*L) ^M^J-^f- 

xc!j(-pi, -p 2 , pi + V2)M{pi)M{p 2 )M{\vi + p 2 |)-P$(pi)-P*(P2) 
25 Ph(k) f d?pid?q T . , T . . „ . , „ . , 

H~-fconst. ) 
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. 1 1 /" d Pi ■ ■ - d Pn L / r, v 

+2/NL2^ („_!)! y C n(Pl>-" .Pn-l.-^n-lJ 

XC^ +1 (- Pl) ■ ■ • , - Pn ,P n )P h ( Pl ) ■ ■ ■ PL(p n -l)M{\P n -i\)M(pn)M(\Pn\) 

x [2P $ (|P n _ 1 |)P $ (p„) + P $ (|P n |)P $ (p„)] 
^ 1 f d 3 Pl ■ ■ ■ d 3 Pn+1 L 

+ 6 -9NlL^T)1 7 ( 27r) 3n + 3 «t(Pl,-,Pn-l,-Pn-l) 

X^ +2 (-Pi, • • • , -p„+i,P„ + i)P L (pi) • • • P L b„-i)X(|P„-i|)X(p„)X(p„ + i)X(|P„ + i|) 
x [3P$(|P„_i|)P$(p n )P*(p„+i) + P*(p n )P*(p n+1 )P4(|P n+1 |)] 
, 25 ^ 2 f d 3 Pl ■ ■ ■ d 3 Pn+1 L 

+ y TNL Z^ („_ i)! y — ( 2 ^)3«+3 — c «(pi>--- ,p„-i,-p„-iJ 

XCn+ 2 (-Pl: ' • • . -Pn+l,Pn+l)PL(Pl) ■ ■ ■ P L (p n -l )-M (|P n -l \)M { Pn )M ( Pn+1 )M (|P n+1 1 ) 

x [p 4 ,(|p Il _ 1 |)p$(p„)P$b„+i) + p$(p„)p$(p„+i)p$(|p„+i|)] 

+6.9NL | ^|^^ G (p 1 ,-p 1 )^ G (-q,q)M(pi) 2 ^( g ) 2 P*(pi) 2 P*( g ) 
^ i f d 3 Pl ---d 3 Pn ^d 3 q L 

+3gNL2_, (n _ 2) , y C n(Pl'--- ,Pn-l,-P„-l) 

xc£(-pi, ■ ■ • , -P«-a, -q, Pn-2 + q)A^(p„_i)>f(|P T[ _i|)M(g)M(|P n _2 + q|) 
xPl(pi) • • • Pl(p»-2) [P*(Pn-i)P«(|P»-i|)P*(?) + P*(p„_i)P$(g)P*(|P„_ 2 + q|)] 

25 ^ 1 f d 3 Pl ---d 3 Pn _ 1 d 3 q L 
+ y TNL ^ 3 (^2)! y (2^ Ctpi,--- ,P»-i,-Pn-i) 

xc£(-pi, • • ■ , -p„_ 2 , -q, Pn-2 + q)A*(p n -i)M(|Pn-i|)A4(g)A't(|Pn-2 + q|) 



xPl(pi)---Pl(p„- 2 ) 



P*(|P„-i|)P*(p„_i)P*(g) 



+P*(p„_i)P*(|P„-i|)P*(|p„-i -q|)+P*(g)P*(|P„-i|)P*(|p„-i -q| 



(34) 



In the above equation, all the higher order correction terms are included in P con st. in large scale limit and P CO nst. gives 
the correction in the bias as A& const . = P CO nst. /Pl(&) oc 1/Pl(/c) oc 1/k. In the pure Gaussian case (/nl = 5nl = 
tnl = 0), the correction in the bias Ab const ./b x a s a function of b\{M) is shown in Fig. [§] From this figure, we 
find that the contributions from the higher order correction terms to the bias parameter are negligible for the objects 
with not so large mass. In case we consider here, the dimensionless power spectrum of the biased objects, which is 
given by k 3 Px{k), is smaller than unity. We could not neglect higher order correction terms any more, when k 3 Px(k) 
becomes order of unity, in which we consider the smaller scales (larger k) and the higher peak objects (larger b x ). We 
also plot the two-loop contribution which is linearly proportional to the non-linearity parameter /nl in Fig. 1101 In 
this figure, we plot \Ab const ./b x \ with /nl = 40 at z = 1.0. We also fix the scale as k = 0.005 h Mpc -1 . This figure 
also shows that the correction terms linearly proportional to /nl are negligible for the objects with not so large mass. 
Here we only show the results up to the two-loop order contributions. However, the higher order contributions than 
three- loops are expected to be small for the case we consider here, because in Fig. |H]the two-loop correction is smaller 
than the one-loop one. We expect that the contributions from the higher order loops which are linearly proportional 
to the primordial trispectrum are also much smaller, from the fact that the two-loop /NL-correction is smaller than 
the Gaussian higher order loop corrections shown in Fig. [TO] Of course, for the more precise discussion about the 
higher order corrections, we need numerical calculations. 
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FIG. 9: A& cons t./&jf as a function of the Eulerian bias 6i(M) for the pure Gaussian case. We fixed z = 1.0 and k 
0.005 h Mpc -1 . The red circles correspond to the one-loop correction and the blue boxes are the two-loop correction. 
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FIG. 10: |Abconst./&xl as a function of the Eulerian bias fei(M) with / NL = 40. We fixed z = 1.0 and k = 0.005 h Mpc 



Finally, we discuss the contributions of the higher order correction terms in the stochasticity parameter. Ref. [l[ 
have discussed that the one-loop term could generate the stochasticity even in the pure Gaussian case. As we have 
mentioned in the previous section, the stochasticity parameter would be a powerful tool to test the relation between 
/nl and tnl and hence the higher order corrections might become "noises" for the test even on large scales. Based 
on the expression for the power spectrum of the biased objects given by Eq. which include the higher order 

corrections, we obtain the stochasticity parameter as 



r(k) - 1 



25 
1 1 



TNL 4/ ^ L ' 6 1 (fc) 2 X(fc) 2 



° <$(p,-p)PL(p) 



(27T) 



2b 1 (k) 2 P L (k) 



d " P ^(p,-p) 2 P L b) 2 



(27T) 



1 ^c£(pi,P2, -Pi - P2) 2 Pl(pi)Pl(p2)Pl(\Pi + PS I 



(2tt) 3 (27r) e 



2/nl j f d 3 Pl d 3 p 2 l 

-^5S-«sCPi -Pi)^ 



6!(fc) 2 P L (fc) 



C2(Pl>-Pl) c 3(Pl>P2:-Pl - P2) 
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xPM^H-Mip^Mdp! + p a |) [2PM + P*(| Pl + pa| 



b 

6.9NL 



(2tt) 3 



(35) 



6 1 (fc) 2 P L (fc) 7 (2tt)3 

Here, we consider the contributions up to the two-loop order. In Fig. 111! we plot the contribution of each term to 
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FIG. 11: Contribution of each term to the stochasticity parameter as a function of the wavenumber, k, with fixing 2 = 1 
and M = 5 x 10 13 ft M©. The red thick line shows the term which is related with the inequality between the non-linearity 
parameters /nl and tnl for /nl = 40 and tnl = 5 x 36/ nl /25. The black solid and dotted lines are respectively for the 
one-loop and two-loop contributions which appear even in pure Gaussian case. The red dashed and dotted-dashed lines are for 
the terms linearly proportional to /nl, which are related with only C2 (negative sign term) and C3, respectively.. 

the stochasticity parameter as a function of the wavenumber k, with fixing z = 1 and M = 5 X 10 13 ft _1 Af Q . The red 
thick line shows the term which is related with the inequality between the non-linearity parameters /nl and tnl for 
/nl = 40 and tnl = 5 x 36/^/25. The black solid and dotted lines are respectively for the one-loop and two-loop 
contributions which appear even in pure Gaussian case. The red dashed and dotted-dashed lines are for the terms 
linearly proportional to /nl, which are related with only (negative sign term) and C3, respectively. For these 
parameters, we find that the relation between /nl and tnl would be obscrvationally checked by large scale survey 
with k < O(10~ 2 ) ftMpc -1 . However, as I mentioned before, the stochasticity includes various uncertainties due to 
other astrophysical processes. The above discussion even including the higher order contributions is still ideal and 
hence the detection of the effect of tnl needs more careful investigation. 



V. SUMMARY AND DISCUSSION 



In this paper, we have derived an accurate formula for the bias parameter with the primordial non-Gaussianity 
parameterized not only by the non-linearity parameter, /nl, but also by <7nl and tnl, by making use of the iPT. The 
scale-dependency of the bias induced from (?nl is the same as that induced from /nl- In this sense, it is difficult to 
distinguish the effect of (?nl with that of /nl- However, we show that the redshift-dependence of the effect of <7nl is 
different from that of /nl and also to obtain a significant constraint for c/nl the higher redshift observations would be 
invaluable. Recently, as an observation of the biased object at high redshift, an ionized fraction in the reionization era 
( 6 < z < 20) through the 21cm observation has been proposed. In Ref. [22|,[23||, the authors shown that through such 
kind of the future observations we could obtain a tight constraint for /nl- We expect that such future observation 
can be also expected to give a tighter constraint for iJnl- 

For tnl , there exist a special inequality between /nl and tnl ■ We consider the possibility of studying this inequality 
through the LSS observations and find that the stochasticity parameter deviates from unity in case tnl > 36/§ L /25. 
We expect that the future wide field survey would be a powerful tool to obtain a constraint for tnl- 
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In deriving the formula with the primordial non-Gaussianity parameterized by <7nl and tnl, we consider the two- loop 
order corrections in the context of the iPT. In usual, these corrections can be considered as higher order effects and 
neglected. Of course, even in case with the pure Gaussian primordial fluctuation there exist higher order contributions 
which may generate the extra scale-dependency of the bias parameter. We investigated these higher order contributions 
in the pure Gaussian case and also in the non-Gaussian case up to the linear order in the non-linearity parameters. 
We find that these higher order contributions can be negligible for the range where the normalized power spectrum 
of the biased objects is smaller than unity, namely, k 3 Px(k) < 1. 

Here, we just analytically derived the bias parameter with the higher order primordial non-Gaussianity. As a 
future issue, hence, we have to check the validity of our formalism with performing numerical simulations and it 
should be interesting to consider the expected constraints for ^nl and tnl in the future surveys. 



NOTE; During the time that we were preparing this manuscript, Ref. 24J appeared on the arXiv. In Ref. 24[, the 
authors focused on the stochasticity due to the primordial non-Gaussianity by using the barrier crossing formalism 
and the peak-background split method. They derived the stochasticity coefficient which is consistent with our result 
in the case with Press-Schechter mass function. 
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